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Abstract 

The study of two-pseudoscalar and two-photon decays for the scalar meson 
nonet below 1 GeV is performed within an effective approach in which the 
scalar resonances are described as (Jaffe's) tetraquark states. The dominant 
(fall apart decay) and the subdominant (one transverse gluon as intermedi- 
ate state) decay amplitudes are systematically taken into account. The latter 
improves the agreement with the experimental data. Possible scenarios con- 
cerning the scalar-isoscalar mixing are discussed. 



1 Introduction 



The interpretation of the mesonic scalar states below 2 GeV is not yet univocal estab- 
lished [TJ 121 El • According to the most popular scenario, one interprets the isovector 
and isotriplet resonances ao(1400) and K(1430) as the ground-state quark- ant iquark 
bound states. The three isoscalar resonances / (1300), / (1500) and / (1710) are 
a mixture of two isoscalar quarkonia and bare glueball configurations (we refer to 
[HSGHniiniElIElin] and Refs. therein). As a consequence, the scalar states below 
1 GeV (a, k, / (980) and a (980)) must be something else, like (loosely bound) 
mesonic molecular states [TUl HI] or Jaffe's tetraquark states [H H2J CS1 El • In this 
work we explore by means of an effective approach the phenomenological implica- 
tions of Jaffe's states, whose building blocks are a diquark (g 2 ) and an antidiquark 
(q 2 ); calculations based on one-gluon exchange support a strong attraction among 
two quarks in a color antitriplet (3c), a flavor antitriplet (3p) and spinless config- 
uration PIE] (color and flavor triplets are realized for an antidiquark). Naively 
speaking, a diquark 'behaves like an antiquark' from a flavor (and color) point of 
view [13J: 

[u,d\*-*s, [u,s]<-+d, [d, s] <— > u, (1) 

therefore out of a diquark and an antidiquark one can build a full scalar nonet, 
whose most appealing property is the reversed mass-order, thus explaining the (al- 
most) degeneracy of the isoscalar state / (980) (whose dominant contribution is the 
tetraquark structure il ss(uu + dd) ,y ) and the isovector a (980) (the neutral one inter- 
preted as u ~ss(uu — dd)"); the a (dominantly ll uudd n ) is then the lightest, in between 
one expects the kaonic state fc(800) {k + interpreted as u ddsu"), which is omitted in 
the compilation of PDG [1*5], but listed in many recent theoretical and experimental 
works ( H51 H21 HBl HS| EDI Ell E21 E3| and. Refs. therein) . 

Evidence for a broad a state is also found, together with a (980) and / (980), in 
the theoretical work of [2~~JI2I~] where a unitarized ChPT is used; as shown in [21*1 "21"* 
a broad scalar state k exists as well. The fact that a full scalar nonet is generated 
within the same approach points to a similar inner structure of the low-lying scalar 
states. Furthermore, the large- N c behavior of the light scalar states indicates large 
non-gg amounts in their spectroscopic wave function |2"""|2""". Support for the Jaffe's 
states has also been found in Lattice calculations (2011301, where the diquark q 2 and 
the antidiquark q 2 are connected by a flux tube. 

In pTT) 13*2"] the scalar states below 1 GeV have also a dominant tetraquark struc- 
ture; mixing among light tetraquark states and heavy quarkonia states, generating 
two scalar nonet of mixed states below and above 1 GeV is described. 

In the present work we intend to analyze the two-pseudoscalar decays of the 
scalar states below 1 GeV when interpreted as Jaffe's tetraquark states; to this end 
we consider the dominant and the subdominant diagrams in the large- N c expansion 
which describe the transition of a tetraquark scalar state into two pseudoscalars 

[121 IH. 
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The dominant diagram, depicted in Fig. l.a, occurs by switch of a quark be- 
longing to the compact diquark with an antiquark of the anti-diquark, thus gener- 
ating two qq objects, which separate as pseudoscalar mesons: q 2 q 2 — > (qq)(qq) — > 
(qq) + (qq). The two pseudoscalar mesons fall apart from the tetraquark configura- 
tion; this decay mechanism was denoted as OZI-superallowed in Ref. [P2] . 

The subdominant diagram, depicted in Fig. l.b, occurs via an annihilation of 
a quark and an antiquark (into one gluon), with subsequent qq creation and two- 
pseudoscalar decay: q 2 q 2 — > qgq — > (qq) + (qq). Although suppressed by a factor 
N c , the fact that the annihilation of the quark- antiquark pair can occur with only 
one gluon as intermediate state (as already noted in Ref. [T2*]). i.e. to order a s only, 
may indeed indicate that the corresponding amplitudes are not negligible. 

In Ref. jT2j the subdominant coupling has not been considered in the decay rates; 
in ^3] it has been introduced as the last step of the analysis in order to improve the 
results of the superallowed decays. In the present work we intend to systematically 
write down the expressions for all scalar-to-pseudoscalars transition amplitudes as 
functions of the strengths of the dominant ("fall apart") and subdominant ("one 
intermediate gluon" ) diagrams, denoted as c\ and C2 respectively. Then, by a fitting 
procedure to experimental known branching ratios, we determine the quantity c%jc\, 
which measures the intensity of the subdominant decay mechanism with respect to 
the dominant one. Two possible solutions with a non-zero value of the ratio c^jcx 
are discussed, which in turn correspond to different isoscalar mixing configurations. 

An interesting feature of the subdominant decay is the strong enhancement of 
the coupling g\_^ K with respect to 9^ o ^ K , thus possibly solving the problem of 
Jaffe's model mentioned in j2I] at the leading (OZI-superallowed, Fig. l.a) order for 
which q 2 _ /o 2 _ = 1 (if no isoscalar mixing occurs, even smaller if mixing is 

introduced), in clear contrast with the result 9\_^ K I 9 % aQ _^ K — 2. 15 ±0.40 reported 
in the analysis of Refs. [2IJI221- The subdominant decay mechanism of Fig. l.b can 
explain the experimental value without introducing explicit KK clouds dressing the 
scalar resonances. 

We then turn our attention to the decay of scalar states into two photons, which 
in line of the previous discussion is described by two decay mechanisms; they are 
analogous to the strong decay diagrams of Fig. 1, where one replaces the two 
pseudoscalar mesons in the final state with two photons. We present the theoretical 
ratios and the phenomenological discussion. 

The paper is organized as follows: in section 2 we write down the Lagrangian 
for the description of the scalar mesons below 1 GeV as tetraquark states; the 
mass term and the two-body strong decays are presented. In section 3 we perform 
a phenomenological study with the available experimental data. In section 4 we 
describe the two-photon transitions and in section 5 we drive our conclusions. 
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2 The model 



2.1 The Lagrangian 

The basic terms are the pseudoscalar fields, collected in the matrix V = 4= Y^h=o B l K (the 
Aj are Gell-Mann matrices), and low- lying scalar fields, collected in the matrix S 
defined as (see Appendix A): 



/ 0\b 



k° 



k+ 



HI 



B 



\ 



J 



h[u,d][u,d] ^[u,d][u,s] d][d, s] 
l[u,s][u,d\ 



§[tt, s][u, s] l[u,s){d,s] 
II 



\ [d, s] [u, d] \ [d, s] [u, s] \ [d, s] [d, s 



(2) 



where in the second matrix the diquark-antidiquark decomposition has been made 
explicit. The states a B = d]\u, d] and f B = i^([u, s][u,s) + [d, s)[d, s]) refer 
to bare (unmixed) states. A mixing of these configurations, leading to the physical 
states a and / (980), is possible and considered below. 

The Lagrangian for the scalar-pseudoscalar interaction reads: 



C 



+ C 1 



(3) 



where (...) denotes trace over flavor. Some comments are in order: 

• In the first term the quantity xp = B ■ diag{m u ,md = m u ,m s } encodes 
flavor symmetry violation. It corresponds to the lowest order chiral perturbation 
theory result [HHl (see also [31] and Refs. therein), to which we refer for a careful 
description. The second term Cf^ ix — — ^ (P ) 2 — zpP°P 8 takes into account the 
enhanced flavor-singlet mass ({7a(1) anomaly) and the octet-singlet mixing, leading 
to the physical states (we follow the notations of [Oj) 



n 



P s cos Bp - P° sin 9 P , rj' = P s sinB P + P° cosB P , 



(4) 



where Bp is the pseudoscalar mixing angle. According to the the standard proce- 



dure [S3 ESI IHZ] we diagonalize the corresponding 77°-?7 8 

139.57 MeV, M K = 



547.75 MeV and Mj 



mass matrix to obtain the 
= 493.677 MeV (the physi- 
957.78 MeV the 



masses of r/ and rj' '. By using M % 
cal charged pion and kaon masses), M v 
mixing angle is determined as Bp = —9.95°, which corresponds to the tree-level result 
(see details in Ref. [HI!)- Correspondingly one finds Mpo = a/ (M| + SM^) /3 + 7p = 
948.10 MeV and z P = -0.105 GeV 2 

• The third and the fourth terms of Eq. (JBJ) refer to the quadratic part of the 
scalar tetraquark nonet; xs — diag{a, f3, f3} and a < (3 takes into account that the 
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diquark [u, d] is lighter than the partners (due to m u = < m s ). In this way the 
masses are given as 

M 2 ao = M 2 fB = 2(3, Ml = (« + /?), M 2 aB = 2a, (5) 

corresponding to the reversed mass ordering. The fourth term C^ ix describes the 
mixing of Ob and details about it are presented in the next subsection. 

It is interesting to note that the same (inverted) mass ordering for the light 
scalars can be obtained in the framework of the 577(3) linear a model j3H] with 
LOi(l)-anomaly. 

It is important to stress that in this work we consider only tree-level expression for 
the decays (see Sections 2.3 and 3) and no mass-shift via loop diagrams is calculated; 
the Lagrangian- masses for the state a (980) and k reported in Eq. © (derived 
from the Lagrangian (J3j)) are therefore the physical masses in our phenomenological 
approach. Similarly, Mf B and M UB of Eq. (jSJ) would be physical masses if no isoscalar 
mixing among as and fs would occur. The physical masses for the states a and 
/o(980) are determined in the next subsection when isoscalar mixing is considered. 

• The last two terms of Eq. (jSJ) correspond to the large- N c dominant (propor- 
tional to ci) and subdominant (proportional to c<i) two-pseudoscalar decays of the 
low-lying scalars when interpreted as tetraquark states. The matrices A % entering 
in the expression (J3J) are the three antisymmetric real 3x3 matrices: 

A 1 = i\ 2 , A 2 = iX 5 , A 3 = i\ 7 . (6) 

Both terms are S77(3)-flavour invariant; the details about the transformation 
properties are reported in Appendix A, where also the other flavour invariant terms 
are listed (Eq. (gUJ)). 

We do not consider in the present work decay terms which break flavour symme- 
try in virtue of a ^ (3 (Eq. (|5|)). Generally flavour-breaking corrections to the decay 
amplitudes are not large and do not change the qualitative picture; furthermore the 
consideration of such terms would imply a too large number of parameters for the 
three-level decay amplitudes. For these reasons the consideration of such terms is 
beyond the goal of the present paper, but it represents a possible future development 
of our study. 

• As described in Appendix A, the S — > PP interaction term c^Sy (A l A^V) (V) 
is also suppressed of a factor N c with respect to the dominant OZI-superallowed one. 
It is coupled to the flavour-blind pseudoscalar configuration (V) = V3P°, connected 
to the physical fields rj and r{ in equation (J3J). It can however occur with at least two 
intermediate transverse gluons attached to P°, i.e. to order a 2 . It's contribution is 
then believed to be smaller than the diagram of Fig. l.b, which takes place at order 
a s . A gluonic amount in the wave functions of the r\ and r( states would enhance 
this channel, but such an eventuality seems not to occur jSHHUIlE]- Furthermore, 
in the analysis of [42J, where the two-pseudoscalar decays of the experimentally 
well-known tensor mesons are evaluated, a good description of data is obtained 
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without an enhanced flavor-blind channel in the pseudoscalar mesonic sector. These 
arguments lead us not to consider the term c^Sij (A l A^V) (V) in the present work. 
According to our view, the further systematical inclusion of this term, which affects 
the couplings involving the rj and rf mesons only, would be then necessary when the 
experimental knowledge on the light scalar meson sector becomes more exhaustive. 



2.2 Mixing 

We discuss the term C mix of Eq. Q which generates mixing between a B and In 
line with the pseudoscalar sector we consider the flavor singlet and octet four-quark 
configurations S = ^/2/3/b + a/1/3<7b and S$ = y/lj^fs — a/2/3<7b and define 

£> S mix = -^"VsSq - z s S S 8 , (7) 

where a octect-singlet mixing and mass-modification for the flavor blind state are 
taken into account; we also refer to Jl] for this point. 

By using Eqs. (j3J) and ((7|) the quadratic part referring to the isoscalar states f B 
and a B and their relative mixing reads: 

C-isoscalar = ~^ ( ff B ^ B ^ ( /b j ' ^ 

where 

(note that M% B and Mj B receive extra-contributions from (J7|) and therefore modified 
from Eq. ©). 

The physical states a and f = / (980) are then given by 



° 1 = ( cos ^s) sm(^s) \ a B 
/o(980) J I -sin(^) cos(^) J \ fs 



(10) 



where 




BVlB 1 = diag{M*, M 2 f J. (11) 



The experimental fact that Mf ~ M ao suggests a small scalar mixing angle 9s- 
Unfortunately, to determine 9s from the masses can be misleading because of the 
large widths of k and a. Furthermore, small variation of the ao and fo masses and 
mixing can lead to different scenarios. We therefore prefer to determine the mixing 
angle form the decay amplitudes (see section 4). 

Nevertheless some interesting considerations can be done; exploiting the relation 
(Q) = Ml + Mj o (from Eq. (jllj) ; the matrix B is orthogonal) one finds 75 = 
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Mj + Ml — 2 ■ M|. The sign of 75 strongly depends on the values of M a and M^. 
One has then to tune the parameter z$ in order to generate the experimental value 
Mf = M ao ~ 0.98 GeV, then the mixing angle 9s is also fixed 1 . We then distinguish 
among 2 possibilities: 

a) 75 > — ► zs < and 9 s < 0; in the upper state fo the bare components ob 
and fs are in phase, while in the lower state a they are out of phase. 

b) 75 < — > 25 > and > 0; the phases are reversed. 

The sign of 9s is important because of destructive or constructive interference 
phenomena in the decays of a and /q. However, as stressed above, it cannot be 
determined from the knowledge of the masses because of the large uncertainties on 
their values. We consider two examples to elucidate this point. Fixing M k = 0.8 
GeV, we have 7s > (and 9s < 0) for M a > 0.56 GeV; vice versa for M a < 0.56 one 
has 7s < (and 9s > 0). Small changes of Mk generate very different result; fixing 
M k = 0.72 (the pole of E791 19 ) one finds 75 > (and 9 S < 0) for M a > 0.277 
GeV ( therefore, if Mk ~ 0.7 GeV the option 9s < seems favoured). We notice 
that the present results differ from those of |14j . 

In Section 4 we will determine the values of the scalar mixing angle 9s from the 
phenomenology and discuss the implications for the masses. 



2.3 Two-pseudoscalar decay amplitudes 

In the following we report the results for the two-pseudoscalar decay rates. We use 
the following notation: for a given decay mode S — > P\Pi (where S refers to a scalar 
state and Pi{P<i) to a pseudoscalar state) the decay width (for kinematically allowed 
decays) is written as 

PS^ Pl P 2 2 ( 

Ts ^ PlP2 ~ ~Mdf 9s ^ p ^ (12) 

where ps^p 1 p 2 is the three-momentum of (one of) the outgoing particle(s) (in the 
rest frame of S): 



Ps^ Pl p 2 = ^f s \/ M s + (MK ~ M%Y ~ 2{Mk + Ml 2 )Ml (13) 

Notice that g s ^ Pl p 2 already includes the charge multiplicities and the symmetry 
factors (see Appendix B). 

In general, but especially when M s < M Pl + Mp 2 , the expression (JT2J) has to be 
modified taking into account the finite width of the resonance, thus integrating over 
the mass of the resonance by employing a suitable mass distribution (like a Breit- 
Wigner one). However, the expressions for the coupling constant are left invariant 

1 The parameter zs can be found from Detffi] = M%-Mj g (see Eq. Hilt ). Indeed, this is a second 
order equation in z$, so strictly speaking two solutions are possible. The solution corresponding 
to the smallest absolute value \0s\ is the one that we take into account. The disregarded solution 
typically induces a very large mixing, in disagreement with the phenomenology (see Section 4). 
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by these operations. We refer to Appendix B for a brief recall on the connection 
between Eq. (fT2j) and the Lagrangian (jSJ). 

In Table 1 we list the coupling constants for the isovector a (980) and the isodou- 
blet(s) k to two pseudoscalar mesons as function of the dominant and subdominant 
decay strengths c\ and c 2 . The expressions in {...} represent the invariant ampli- 
tudes, eventually multiplied by a factor 2 for decays into identical particles, while 
the coefficients in front of the parenthesis account for charge multiplicities and sym- 
metry factors (see Appendix B). 

In Table 2 we report the results for the bare states <7b and although not 
physical when S ^ 0, the expressions are easier to read and allow to understand 
the role of the subleading decay mechanism of Fig. l.b. It is important to note that 
the decay mode fp — > KK is significantly enhanced with respect to ao — > KK when 
> 0. This fact shows that a non-negligible and positive ratio c 2 /ci can explain 
why the /o — > KK coupling is larger than the ao — > KK one. This point will be 
discussed in the next Section. 

When considering the mixed physical states a and /o(980) defined in Eq. (JTUJl 
the coupling constants are modified as follows: 

g*-,p 1 p 2 = {g* B ^PiP 2 ■ cos(B s ) + g fB ^ PlP2 ■ sm(6 s )} 
£/o-PiP 2 = {-9*b-+i\p» ■ sin (#s) + 9f B -*PiP2 ■ cos(6 s )} (14) 

Because of their relevance, we report in Table 3 the explicit expressions for the 
decays of a and /o(980) into titi and KK as derived from Table 2 and Eqs. (jUj) . 



Table 1. Decay coupling constants for a and k 



S -> PxP 2 


9s^PiP 2 


a — > 1(K 




V2c 1 + j= 2 c 2 ) 




a -> TTi] 


VT- 




v2 cos(6 l p) + sin 


(9 P )] + ^ 


cos(0 P ) - V / 2sin(0 P )] | 




a — * 7CT]' 




{-> 


cos(6» P ) - v / 2sin(6'p)] + y§c 2 [v^cos^p) + sin(0 P )] J 




k — > ttK 




L V2c 1+ ^ C2 ; 


• 


k — > Kr] 


vT-{ 


~V 3 1 


cos(0 P ) - v^sin^p)] - <J\cz [cos(# P ) + 2 v / 2sin(#p)" 


} 


k -> Krf 


Vi-{ 


f^r 
~V 3 1 


V2cos(0p) + sin(0 P )] + ^±c 2 [2\/2cos(#p) - sin(# P )" 


} 
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Table 2. Decay coupling constants of the bare states a B and /p 



s — > 


PtPo 
1 l 1 2 


9S^PiP 2 


u B 


7T7T 


v /|-{2c 1 + 2c 2 } 


(J B -) 


KK 


• {c 2 } 


ctb - 


-> rjr] 


\Zl • {f(- c i + c 2) ( cos (M - V2sin(^ P )) 2, 




b 


'/ '/ 


^1 ■ {f(-ci + c 2 ) (V2cos(0 P ) + sin(0p)) 2 ^ 




ob - 




y/l- {|(-ci + c 2 ) (2 v / 2cos(20p) -sin(2^p)) 


} 


!b~ 


* 7T7T 


a/|-{V2c 2 } 


fs } 




V2-{(^c 1 + ^c 2 )} 


f B - 


-> 77?7 


^1 • { | (9V2c 2 + y/2{8d + c 2 ) cos(2fl P ) + 4(- Cl + c 2 ) 


sin(29 P ))} 


fB~ 


► r/'r/' 


■ {| (9^2 - V2(8 Cl + c 2 ) cos(20p) + 4( Cl - c 2 ) 


sin(29 P ))} 


f B - 




V 7 ! • {i (4( Cl - c 2 ) cos(20p) + v/2(8 Cl + c 2 ) sin(2^ P )} 



Table 3. Decay coupling constants of the physical states a and fo 



S -> PxP 2 


9s^PiP 2 


0" — > vTvT 


[3 

V 2 


{2 Cl cos(0 s ) + C2(2cos(0 s ) + \/2sin(0 s ))} 




x/2- 


{ ^2 Cl sin(fl 5 ) + c 2 (cos(0 s ) + ^ sin(0 s )) } 




/o -» kit 


>/§•{ 


-2ci sin(0 s ) +c 2 (-2sin(0 5 ) + v/2cos(0 s ))} 


h^KK 


v/2-{ 


v/2 Cl cos(0 s ) + c 2 (- sin(0 s ) + ^ cos(0 5 )) 





3 Strong decays: numerical results 

The data about the decay widths reported in are still not complete. The reso- 
nance a is very broad (M a = 0.4-1.2 GeV, T = 0.6-1 GeV), k is still not yet listed (in 
the recent result of (201 the pole is found at 0.76 GeV, the width is very large). The 
masses for the states a (980) and / (980) are well known, M ao = 948.7 ± 1.2 MeV 
M fo = 980 ± 10 MeV, but the widths are not, between 40 and 100 MeV for both 
states. The presence of near-threshold KK decays complicates the experimental 
and theoretical analysis (see discussion in and below). 
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Information about (some of) the coupling constants g 2 s ^ Pl p 2 has been extracted 
directly from experiment [201 EH 1221 ESI- in the case of a (980), for instance, the 
quantities g 2 Q ^ v , 9 2 Q ^ K and M ao are free parameters of the meson- meson scat- 
tering amplitudes (generally the Flatte distribution is used jlSHS]). They can be 
determined by fitting the theoretical amplitudes to the experimental ones. How- 
ever, as explained in [IH], relatively large differences in the absolute values of the 



coupling constants gl ^ nv and found in the literature. Fortunately, the 

ratio o?„ ._,/ Q 2 shows a stable behavior. 

Furthermore, the presence of strong coupling constant for the near-threshold 
decay mode g ao _,~K K affects also the channel a — > tct] rendering the experimental 
width narrower than the theoretical result obtained using Eq. (|12|). see the discussion 
in [25] . 

In the following we will compare our results with the analysis of experimental 
results of Refs. [201 12H 122] 5 where the ratios of coupling constants for various decay 
channels are deduced; it should be stressed that the results of [201 ED 122] are not 
yet conclusive and depend on the choice of the form factors and other assumptions, 
for which we refer to the above cited works for a careful description. 

The ratios of coupling constants for the resonances ao(980) and /o(980) as re- 
ported in [201 [2D [221 are: 



a 2 — a 2 — 2 

f l~* KK = 4.21 ± 0.46, S ^ KK = 2.15 ± 0.40, = 0.75 ± 0.11. (15) 

dfo^wTT 9 ao ^KK Qao^KK 



These three results are indeed (at least qualitatively) common to (almost) all the 
analyses (see [13] and Refs. therein, and the recent experimental analysis of [431 14T] ). 
The largest uncertainty is about the crossed ratio 9 2 o ^ K / 9 2 ao ^K K 'i however, the 

enhanced coupling of the fo — *> KK mode with respect to a — > KK is also a stable 
result. As we will see, all three ratios in (|To'|) are all compatible with a sizable C2/C1. 

Before considering a nonzero value for Cijc\ , we analyze as a first step the case 
C2 = 0, corresponding to the original work of ^2] and described in the recent paper 
of [21]- We fit the only free parameter 0$ to the first and the second branching 
ratios of Eq. (f]3j) ; the ratio g 2 ^ nr) / '5^^ = 0.50 does not depend on S , hence 
not included in the fit (in [21. it is 0.4 because of a slightly different pseudoscalar 
mixing angle; here it is Op = —9.95°, while in [21] the value 6p = —17.29° is used). 

We find the following value for 0s (we refer to this Solution A): 

V 2 

Sol. A: (cjj = 0), S = ±21.6°; (|- = 5.17) (16) 

corresponding to the values 

a 2 — a 2 — 2 

y Jp™L = 4.26, ZJpBL = 0.86, = °- 50 - ( 17 ) 

9f ^TTTv 9 ao ^KK Qao^KK 
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The large and unsatisfactory \ is generated by the mismatch between experiment 
and theory for the ratio g 2 _^ K l g 2 ao ^-^ K - Two scalar mixing angles are reported in 
Eq. (fTtjj) because the fitted quantities are symmetric for 9 S — > —9 S when c 2 = 0. 
Our theoretical ratio for q 2 - F7 „/q 2 f as function of 9s is deduced from Table 3 

JO — >KK ' JO >7f7T " 

and reads for c 2 = 0: 

g JpEH = 2 -cot\6 s ). (18) 

We notice that in [21] the theoretical ratio is <7^ o ^= i< -/<?/ _ >7r7r = |cot 2 (^ 5 ); for this 
reason the mixing angle found in j2I] is 9$ = ±15.9° and differs from our result. The 
extra-factor 1/2 present in the result of |2I] could not be verified when evaluating 
the traces of Eq. Q . 

We now consider a non-zero c 2 : in principle one could determine the ratio 
c 2 /ci from the quantity 9a ^nr 1 /9a ^'KK a ^ one - m ^ ac ^; the theoretical expres- 
sion depends only on c 2 /ci. When c 2 = the result as derived from Table 1 
is 0.5, then gl o ^ nr] / g^^ R increases very slowly for increasing c 2 /ci. The value 

02o-W 'dl^KK = °- 6 is reached when c 2/ci = 0.6, but the value g 2 ao ^ nv /g 2 ao ^ K = 
0.75 corresponds to c 2 /ci = 2.35. Because of this strong sensibility it is not practi- 
cable to determine c 2 /ci in this way. However, a sizable c 2 improves the agreement 
with the experiment. 

The x 2 -method is then used to find the free quantities c 2 /ci and 9s which cor- 
respond to the best description of (JT3J) (Solution B): 

Sol. B: c 2 /ci = 0.62, 9 S = -12.8°; (^- = 0.65 < 1). (19) 



The theoretical ratios evaluated with the parameters of Sol. B are: 
a 2 — a 2 — 2 

Jpl££ = 4.21, ZJp^]L = 2.28, = 0.60, (20) 

9f ^T„T ®a^~KK ^ao^KK 



thus in good agreement with the results of (fTKj) . The inclusion of the subdominant 
decay diagram of Fig. l.b leads to a clear improvement of all three ratios without 
adding KK contributions to the wave functions of the resonances. Notice that 
within this Solution a negative value for the scalar mixing angle is preferred. 

In [2U 122] a large ratio g 2 a ^ K / g 2 ^,^ = 0.6 ± 0.1 is deduced from analysis of 
— * 77r°7r° (such a large ratio can explain a problem in reproducing the overall 
normalizations; other explanations are however possible ^HIIM])- With the param- 
eters of Sol. B (fTTJj) one finds a very small ratio: g 2 _^ K I g 2 ^ n = 4.8 • 10~ 7 . Such a 
small value is caused by the destructive interference between the <jb and Jb compo- 
nents (see Table 3). If the value of (22] should be confirmed by future experimental 
analyses, our Sol. B should be rejected. 
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We make a second fit by adding to the data in (fT5|) the value g 2 a ^ K l gl^ nw = 
0.6 ± 0.1. The minimum of x 2 is found for (Solution C): 

Sol. C: c 2 /ci = 0.89, 9 S = 35.8°; (y = 2.04). (21) 

The results derived from the Solutions A, B and C are summarized in Table 4, 
where also other coupling ratios are presented. In Sol. C a large branching ratio 
KK/nn for the a resonance is found (0.65), however g 2 ^— K /g 2 gets worse 

(1.12, generating a large \ 2 'i it is however still larger than 1). 



Table 4. Comparison of coupling ratios with the analysis of experimental results 

of Refs. El 122 2 



Ratios of g 2 


Sol. a (HHD 


Sol. B dm 


Sol. c (EH) 


Analysis of Refs. [21, 22 


q 2 lq 2 - 


0.50 


0.60 


0.63 


0.75 ±0.11 


n 2 In 2 

y f {) ^KK' JO *7T7T 


4.26 


4.21 


4.35 


4.21 ±0.16 


q 2 - lq 2 - 


0.86 


2.28 


1.11 


2.15 ±0.4 


n 2 In 2 


1.01 


0.16 


0.05 




9fo-*rmi 9f -nnr 


1.78 (1.57) 


1.35 


2.30 


< 0.33 


q 2 - lq 2 


0.10 


4.8 • 10" 7 


0.65 


0.6 ±0.1 


q 2 lq 2 


0.30 (0.02) 


0.05 


0.10 


0.20 ±0.04 


dk^TTK/ 9 (J — >7T7T 


1.16 


0.78 


0.58 


(2.14 ±0.28) to (1.35 ±0.10) 




0.17 


0.12 


0.11 


0.06 ±0.02 


9k->ri'K/ 9k->nK 


0.16 


0.006 


~ 


0.29 ±0.29 



Some comments are in order: 

a) The first three ratios refer to a (980) and to /o(980), which have been studied 
in details in the literature for what concern the KK, irrj and tctc channels. For these 
reasons our preferred solution is Sol. B (|19jh which is in good agreement with these 
three experimental ratios. 

The other ratios refer to the broad a and k states or to channels far from threshold 
for a and /o, therefore are not free of ambiguities. 

b) The ratio g 2 -j( K / 9%-^ ls ^ ne marn difference between Solutions B and C, 
which also generates different values for the scalar angle, negative in the first case 
and positive in the second. We notice also a large discrepancy of all three Solutions 
for the branching ratio gj ^ m / gj ^ nw , whose experimental value reported in [21] is 
smaller than the theoretical results for Jaffe's tetraquark states. Future checks on 
these two quantities may help to disentangle the nature of the light scalar mesons. 

2 The results for Sol. A are equal for the two mixing angles 6s — ±21.6°, with the exception of 
two ratios: g'j g _ rjri / gf ^ n7T an d 9 2 G ^ nn l '9%-nnn f° r which the first result corresponds to 9s = —21.6°, 
while the secon (in parenthesys) to 9 s = 21.6°. 
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c) In Solutions B and C the ratio C2/C1 is not small. The inclusion of the 
decay mechanism of Fig. l.b (which can occur with one gluon as intermediate 
state) represents an improvement for the phenomenology of Jaffe's tetraquark states. 
Indeed, this result has been anticipated in [J4|, where the quantity Cijc\ varies 
between 0.7 and 1 (for comparison, the parameters a and b in ^3] are 2ci and — C2 
in the present work). 

In particular, one can notice that the first three entries of Table 4 are improved 
in Sol. B when compared to Sol. A (corresponding to c-i — 0), thus a non-zero and 
sizable ratio c<ijc\ is favoured within the here presented phenomenological analysis. 

In Sol. B the ratio C2/C1 = 0.62 lies in between 1/N C = 1/3 and 1. Although 
larger than 1/3, it is still clearly smaller than 1. Sol. B is still in agreement with 
expectations from large N c considerations. In Sol. C one has Cijc\ = 0.89, which is 
of order 1. Such a value would imply a substantial violation from large N c ; this fact 
constitutes a further hint in favor of Sol. B. 

d) The coupling g k ^ K is smaller then gl_, n7T in solution B and C ; this is not in 
accord with the experimental result describing such a wide k resonance (curiously, 
Solution A is better in this respect). Two experimental values are reported in Table 
4 (see PlJ), the second one does not represent such a large mismatch with the 
theoretical values. It is noticeable that a similar problem exists when describing the 
scalar sector above 1 GeV (|9, and Refs. therein): the experimental large width of 
the state i^o(1430) cannot be explained theoretically. 

e) The knowledge of the mixing angle does not allow us to deduce M a and M k . 
However, we can consider reasonable values for M a and determine M k for which 
the mixing angles of Sol. B and C are realized: if we take M a = 0.45 GeV, then 
a scalar mixing angle 9s = —12.8° (Sol. B (fT9"j) ) corresponds to M k = 0.69 GeV, 
while a mixing angle S = 35.8° (Sol. C (JUJ)) to M k = 0.944 GeV. Similarly, if 
M a = 0.55 GeV the corresponding couples of values (9 S = — 12.8°, M fe = 0.735 
GeV) and (9 S = 31.8°, M k = 0.953 GeV) are found. The present values for the k 
pole favor a light mass between 0.7 and 0.8 GeV [P9*l |2*U] , therefore Sol. B is in 
better agreement with the data. (However, the analysis of j3H| points to a slightly 
heavier kaonic state: M k = 841±30^73, therefore caution is still needed when driving 
conclusions) . 

/) Keeping in mind the remarks at the beginning of this section, we nevertheless 
discuss some full widths as (naively) calculated from Eq. (|T2|) ; in (33] various values 
for gl ^ nr) are reported (= g ■ 8irM% in [32] )• The results vary between 5 and 10 
GeV 2 . We take for simplicity a value in between: fl'^ _ >7r _ = 7.5 GeV 2 . For Sol. B such 
a value corresponds to C\ = 1.32 GeV, resulting in (values in MeV) r^^^ = 136, 
r ao ^7rr; = 98, I" ' a ^im = 795 and Yk^K-w = 251. For Sol. C we have c\ = 1.19 GeV, 
implying T fo ^ n7T = 58, T ao ^ v = 98, T a ^ n7T = 1013 and T k ^ Kn = 241. Clearly 
the widths change substantially when varying g 2 ^^,^ in the mentioned range. The 
qualitative picture emerging is the presence in both cases of two broad states a and 
k with r^Tt-Tr > T^^,,- (but the latter not so broad as desired), and two narrower 
partial decay widths for f and a ; notice that the ordering Tf^^ > T ao ^ nr] in Sol. 
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B is reversed in Sol. C. The inclusion of a mass distribution for the calculation of 
the decay widths, thus also evaluating the KK modes for a and /o, is planned as 
a future step but beyond the goal of the present work. 

g) Before moving to the two-photon decays, we wish to remind what has been 
considered and what omitted in the performed study on strong decays. The two 
interaction terms of the Lagrangian (J3j) are ST/^ 3) -invariant and correspond to the 
diagrams of Fig. 1; other flavour symmetric terms, listed in Eq. (|4(J|) . and direct 
S'[/y(3)-breaking interactions have been omitted. They are supposed to represent 
corrections to the here presented scenario(s) (large ^-suppressed terms and flavor 
symmetry breaking corrections) eventually useful for a more quantitative study. By 
performing a fit leading to Sol. B we used as experimental quantities the values 
reported in [2H1 ED 122]; as already noticed above, these three results, although 
qualitatively similar to other works on light scalars, depend on particular assumption 
of form factors. Furthermore, the errors are only statistical and not systematic. Even 
more caution is needed with the other values reported in the right column of Table 
4. Being aware of these limitations of both theoretical and experimental origin, we 
however intended to focus on a particular and interesting aspect of phenomenology 
of four-quark states encoded in the two diagrams of Fig. 1. 



4 Two-photon decays 

We now turn the attention to the 77-decays of the Jaffe's tetraquark states. As 
for the strong decays we consider two analogous channels, where one has photons 
instead of mesons as final states in Fig 1. The Lagrangian reads 

C em = cpSy {A L QA j Q) F% - cfSii (AAQ 2 ) F 2 V , (22) 

where Q = e ■ diag{2/3, —1/3, —1/3} is the charge matrix (e = y/4ira is the electron 
charge, a ~ 1/137) and F^ v = d^A v — d u A^ the electromagnetic field tensor. Note 
that the convention for the relative sign of the leading and subleading terms is the 
same of Eq. (j3J) (see also Appendix A). 
The decay width into two photons reads 

r w = 16vr« 2 M| • g 2 ^, (23) 

where the nonzero contributions obviously correspond to S = a[j, er, fo- The coupling 
constants for a[j and for the bare states erg and fs are deduced from and read: 

2cJ< + c 77 4c 77 + 5c 77 _ 2c 77 + 7c 77 

When considering the mixed physical states a and /o(980) defined in Eq. (J1U|) the 
coupling constants are modified as (see also Eq. Q14|)): 

# CT -> 77 = 9<x B -^n ' cos (#s0 + fi'/ s ^77 ' sin(0 s ), 
57o-77 = -9*B-~n ' sin(0 s ) + #/ s ^ 77 • cos(6>s). (25) 
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The experimental results for the decay width of ao and fo are given by 

r /o-77 = O-SS-als KeV > r a^ 77 = 0.30 ± 0.10 KeV. (26) 

The experimental value for r a o^ 27 is not reported as an average in however, 
it is in accord with the quoted averages r a o 727 • (T V7T /T to t) = 0.24+q;q7 KeV and 
Fk K /T V7T = 0.183±0.024. Combining these two results one finds T a o^ 77 = 0.28±0.10 
KeV, well compatible with (|26|). In [Slj it was shown that two-photon widths about 
0.3 keV are compatible with a tetraquark nature of the states. 

Expressing the results of (j2T?j) in terms of squared coupling constants we find 



9%_ 



r 



77 _ ^Qq^77 
/o^77 9fo-*TI 



0.77 ±0.48 (27) 



(where M ao = Mf = 0.98 GeV and an average error of 0.115 KeV for r/ _^ 77 have 
been used). 

Large N c results are still valid by invoking vector meson dominance, for which 
the decay into two photons occurs via two virtual vector mesons. Then, we first 
discuss the case c 2 7 = (Solution A), where only the dominant contribution is 
considered. By using the mixing angle 6 s = 21.6° (see Eq. (jl6Jl ) we find the totally 
wrong ratio <7^o_ >Ty /<7/ _> 77 = 724.7, while the other option 6 S = —21.6° implies 
9a -* /#/o->77 = 2-3, which is still unacceptable when compared to the experimental 
result of (|27|). Even by modifying the scalar angle 6 s the situation is not improved; 
one finds a minimum for the ratio q 2 I q 2 f = 1 at 6s — —70.5°. Such a value 
for the scalar mixing angle is however ruled out by the phenomenology of the strong 
decay analyzed in the previous section (a small mixing angle 6 s is common to all 
three analyzed scenarios and is in accord with the mass degeneracy of ao(980) and 
/o(980)). Therefore, also the two-photon decay shows that the dominant decay 
mechanism, corresponding to a switch of q and q analogous to Fig. l.a, is not 
enough to describe the experimental data. 

Let us then consider the case of a non-zero op ■ We keep the scalar mixing angle 
9s as found in Sol. B and C respectively and we determine the ratio cp 1 /cj 1 in order 
to reproduce the experimental result of Eq. (J2TI) 3 

c 77 

Sol. B (Eq. (JED, 9 S = -12.8°) = 0.73 (28) 



77 
2 



Sol. C (Eq. (HQ), 6 S = 35.8°) = -1.04 (29) 



„77 



In both large ratio cj 1 /cj 1 is found; while in Sol. B it is still safely smaller 

than unity, in Sol. C a large (and negative < — 1) ratio is obtained. 

3 Once the ratio is determined, one can fix the strength of of 1 to reproduce r/ ^ 77 (Eq. (|26[l '). 
The other width r/ ^ 77 is then also correctly described. 
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Furthermore, in Solution B we find a noticeable correspondence between the 
electromagnetic and strong transitions: c£ '/of = 0.73 ~ c 2 /ci = 0.62. The two- 
photon decay, which in the framework of vector meson dominance is mediated by two 
virtual vector mesons, differs from the two-pseudoscalar decay only in the final stage, 
therefore the contribution of the subleading diagram with respect to the leading one 
is naively expected to be of the same magnitude (and with the same sign). On the 
contrary in Solution C we have c^/c 77 = —1.04 << c 2 /ci = 0.89 : this fact would 
mean very different contributions of the subleading diagram in the decay amplitude 
into two vector mesons (then converting into two photons) from the corresponding 
transition into two pseudoscalar mesons. The composite approach used in this work 
does not allow for a more microscopic insights to deal rigorously with this issue, 
therefore this discussion does not represent a proof in favour of Sol. B. We simply 
limit to notice that within Sol. B one has c^/c 77 ~ C2/C1, which on the contrary 
does not take place in Sol. C. 

The experimental situation concerning a —>■ 2 / y is less clear; no average or fit is 
presented in ^3] , however two experiments listed in find large 77 decay widths: 
3.8 ±1.5 keV and 5.4 ±2.3 keV, respectively. In a footnote it is then state that these 
values could be assigned to / (1370) (actually, in a older version of Particle Data 
Group these values were assigned to the resonance / (1370)). It is not clear if 
the 77 signal comes from the high mass tail of the broad a state or from /o(1370) 
(or even from both; in such case the experimental result would represent the sum, 
i.e. an upper limit for both resonances). Furthermore, the application of Eq. ()23|) 
is too naive for the broad a state; in fact, due to the third power of the mass of the 
decaying resonance, the two-photon decay width is strongly influenced by a large 
width, especially from the right mass tail (the precise form of a mass distribution 
would be necessary for a more precise and quantitative analysis). Here we simply 
report the results for the ratio 9^0^/9^^ i n the t wo scenarios: g^p / ^-^77 * s 
0.83 in Solution B ([2~%|) and 0.42 in Solution C (J2HJ). However, if a two-photon width 
between 3 and 5 keV should be entirely assigned to the cx-resonance as recently 
discussed in Ref. the interpretation of the a as a tetraquark state would be 

problematic, because in the four-quark scenario the 77 decay width of the a would 
be of the same order of the giq(980) one (~ 0.3 keV) , as the above reported coupling 
ratios g 2 a o^ / 'fl£_> 77 suggest. Future work on the o — > 77 transition is needed both 
theoretically and experimentally in order to clarify this crucial issue about light 
scalar mesons. 

5 Summary and Conclusions 

The Jaffe's four-quark states have appealing characteristics to be good candidates 
for the description of the scalar nonet below 1 GeV; this paper aimed to analyze this 
possibility by studying strong and the electromagnetic decays of scalar resonances. 
We first considered the strong decays of the light scalar mesons below 1 GeV. 
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Beyond the OZI-superallowed decay, in which the scalar tetraquark state falls apart 
into two pseudoscalar mesons as depicted in Fig. l.a, the next to leading order in 
the N c expansion, which occurs with (only) one intermediate gluon as shown in Fig. 
l.b, has been systematically taken into account for all decay coupling constants 
and presented in Tables 1, 2 and 3. The two decay mechanisms are described 
by corresponding terms in an effective composite interaction Lagrangian and are 
parametrized by two strength parameters C\ and c 2 respectively. 

As a first step we studied the case c 2 = 0, in which the diagram of Fig. l.b is 
switched off. The large ratio 9 2 _^ K I 9 2 aa _^ K — 2.15 ± 0.40 > 1 reported in the 
experimental analysis of [22] cannot be described because the theoretical ratio is 
< 1 for C2 = 0. We referred to the case c% = as 'Solution A'. 

We then allowed for a non-zero ratio cz/ci. by fitting c<ijc\ and the scalar- 
isoscalar mixing angle 9s to the three values reported in Eq. ([To]) deduced in the 
analysis of [22; and involving the KK, iri] and tttt channels for the a (980) and 
/o(980) states, a sizable ratio cijc\ = 0.62 is obtained. The subdominant transition 
of Fig. l.b is non negligible in the present analysis and improves the theoretical 
description of the mentioned decay modes of ao(980) and /o(980). The mixing an- 
gle 9s = —12.8° is small and negative: the bare components o B = u uudd" and 
fs = u ss(uu — dd)" are out of phase in the physical state a and in phase in the res- 
onance /o(980). We called this scenario 'Solution B', which represents our preferred 
solution, as discussed in Sections 3 and 4. 

Within Solution B one finds a very small ratio g 2 ^ KK /g 2 ^ 7Tn ~ 10~ 7 , while the 
result reported in the analysis of Ref. [21] reads 0.6 ± 0.1. Solution C has been then 
built by including this ratio when minimizing the x-squared. As a result, agreement 
with 122] is obtained for the branching ratio KK/ittt, but the ratio q 2 , -^^Iq 2 

1 ' ° ' Jo— *KK' ^ao^KK 

is found to be 1.11 within Solution C, which is, although larger than 1, clearly 
worsened with respect to Solution B. Future experimental check on g 2 ^ KK / g 2 ^ n7T 
may help to clarify this point. The mixing angle in Sol. C is positive: 9g = 31.8°, 
i.e. the bare components are in phase for o and out of phase for / (980). The 
sign difference of 9s is responsible of such different results concerning the quantity 

— >7T7r* 

The results of Solutions A, B and C are summarized in Table 4, where also other 
ratios of squared coupling constants are listed. 

We then considered the two-photon decay of the scalar resonances. The diagrams 
for such transitions are similar to those of Fig. l.a and l.b but with photons instead 
of mesons as final states. The corresponding intensities are parametrized by of and 
c 2 7 respectively. 

The discussions follows the same line as before: the case c 2 7 = (Sol. A) is 
analyzed first. The experimental result r a o^ 27 /Tj ^ 2 7 = 9 2 a a^ n l 'g} ->^ = 0.77±0.48 
cannot be obtained if cj 7 = 0, thus again pointing to a non-zero contribution of the 
subleading diagram with one intermediate gluon. 

We then analyzed Solutions B and C: by taking the scalar mixing angle 9s from 
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Solution B or C the ratio cP/cf has been determined in order to reproduce the 
experimental value g*o_^ /fl'/ _ i . TY = 0.77 ± 0.48. In case B one has the noticeable 
correspondence cP/cf ~ c 2 /ci, that is the contribution of the subleading decay is 
similar (and with the same sign) in strong and electromagnetic transitions. This fact 
does not take place in Solution C where the two ratios are large but with opposite 
sign. 

The interpretation of the light scalar states as tetraquark objects is a viable and 
interesting possibility. The inclusion of the subleading decay mechanism of Fig. l.b 
improves the description of experimental data for both strong and electromagnetic 
transitions within our Solution B. 

Future work is however needed: on a experimental side a clear experimental 
establishment of the resonance k and a better understanding of the parameters 
describing the broad states a and k would constitute a decisive improvement. 

We would like to mention possible future theoretical studies: the decays of heav- 
ier states, such as tensor mesons or scalar mesons above 1 GeV, into two light 
scalar mesons as era, kk... should be evaluated by explicitly taking into account the 
four-quark nature of the light scalars. Such a analysis would therefore extend and 
complete the ones in (211121 ■ in particular, the fact that the tensor meson nonet is 
well established and well known experimentally would offer a suitable environment 
for such a study. Furthermore, the two-body decay of excited pseudoscalar states 
(eventually mixed with a pseudoscalar glueball) would also constitute an interest- 
ing development. In the recent work of [53] attention is focused on electromagnetic 
decays of the type S — > Vj, which may also be helpful in disentangling the nature 
of the light scalar states below 1 GeV. 
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A Tetraquark scalar nonet 

The wave function for a "good" diquark is [T3] 

\qq) = |space: L = 0) |spin: S = 0) [color: 3c) |flavor: 3p) . (30) 

The state \qq) is antisymmetric by exchange of the two particles in accord with the 
Pauli principle; the parity quantum number isP=(— 1) L = +1. In the following we 
consider the flavor decomposition for a diquark (and for an antidiquark). To this end 
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we define the vector q such that q f = (u,d,s); the antisymmetric 3j? decomposition 
of a diquark is described by the antisymmetric matrix D: 

D=^j\{qi-q t 2 -q2-q\) = ^A % (31) 

where the matrices A 1 are the three antisymmetric 3x3 real matrices of Eq. (JBJ) 
and the quantities <fi are given by 

<Pi = (uid 2 - u 2 d l ) = [it, d] , (32) 

^ = y| («is 2 - «2Si) = [«, s] , (33) 

^ = \[\ ( rf i s 2 - tfesi) = [d, s] . (34) 
The anti-diquark is then described by the matrix 

Dt = jp t (Ay = -WiA 1 (35) 



with cp 1 = a/1/2 [it, <i] , <£> 2 = a/1/2 [it, s] and 9? 3 = a/1/2 [5, s] . 
The matrix D has the following transformation properties: 

Sf/(3)-flavor: D -> C/ G SU(3), (36) 

P-parity: D -»• D, (37) 
C-Charge Conjugation: P — ► D^; (38) 

The flavor transformation (|3fijl. corresponding to g — > C/g, follows directly from the 
definition (|3T|). 

The matrix 5 defined in Eq. (J2J) and entering in Eq. (J3J) is given by <Sjj = y^j. 
In a generic Lagrangian term involving linear tetraquark scalar states (i.e. one 
matrix D and one matrix D*) one has a schematic form like 

(.../;.../; ...) ^ (...,r...(,F) ...) (39) 

where dots indicate operation with matrices, such as multiplication, traces, etc. The 
form on the l.h.s. is useful to check the iS77(3)-flavor transformations, while the form 
in the r.h.s. allows the calculations of the decay amplitudes (see Tables 1, 2, 3 and 
Appendix B). When describing the scalar — > two-pseudoscalar decay we then have 
the following possible flavor-invariant terms: 

£s^pp = -c 1 (DV t D^V) + c 2 (DD^V 2 )+c 3 (DD^V)(V) 

+c 4 (Drf)(V 2 )+c 5 (Drf)(V}\ (40) 

where V = 4j Y^=o P l ^i- The first two terms corresponds to the two diagrams of 
Fig. 1, i.e. to the last two terms of Eq. Q (by applying the operation (J3~9"j) the 
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equality is easily seen) described throughout all the work. The third term is also 
suppressed by a factor N c only with respect to the OZI-superallowed decay, but 
the transition can occur with at least two transverse gluons as intermediate states 
attached to the pseudoscalar flavour- singlet state. The fourth and the fifth terms 
are further suppressed. 

The Lagrangian £$^pp is clearly charge and parity invariant (we recall that 
V ^UVU^ , V — > — V and V — » V under flavor, parity and charge conjugation trans- 
formations). 

Formally the two-photon Lagrangian (|22|) can be obtained from Eq. ()4(JJ) by 
replacing V — > QF^ U1 Ci — > c] 1 (and then contracting the Lorentz indices). In 
virtue of (Q) = the first, the second and the fourth terms 'survive' the replacing; 
then, neglecting the latter we get the Lagrangian (|2*2*|) (by invoking vector meson 
dominance, large N c arguments can be applied also for the two-photon decays). 

When considering the quadratic term for the scalar-tetraquark states care is 
needed; let us introduce an extra-index a = 1,2 with D a = tp\A % . The elements 
Sij (Eq. ((21)) are formed by the two corresponding diquark (pf and antidiquark Tp a -\ 

= ^\Tp a -. Let us consider the flavor-invariant term (D 2 )^ (d 2 (D 1 )^ ; we 

aim to show that this is nothing else but 4 (S 2 ) (proportional to the mass term of 
Eq. (JHJ) in the flavour-symmetric limit): 

(D 2 y) (d 2 (o i y) = <^> <^<> 

= 4(riy>lM<p 2 j)5ij5 kl = ASuS H = 4 (S 2 ) . (41) 

Then, in Eq. (JSJ) modification from the flavor symmetric limit are parametrized by 
the matrix xs = diag{a, j3, 0} and a < (3. 

In Section 2.2 (Eq. ((7|)) we considered extra-terms in the scalar singlet-octet and 
singlet-singlet channels. The singlet-singlet channel, proportional to S$, is flavour- 
symmetric and corresponds to l^D 1 (Z) 1 )^ (^D 2 (D 2 )^ . The singlet-octet mixing 

introduced in Eq. (|7J) breaks flavour invariance. Here we notice that a flavour 
symmetric term, also including singlet-octet mixing (but affecting also the isovector 
and isodoublet states) can be constructed for the quadratic scalar sector: 

(D 1 (D 1 ) 1 D 2 (D 2 )) = S u S kl (A i A j A k A l ) . (42) 

Further study on the mass sector for tetraquark scalar states including this term 
could be interesting. 

B Strong coupling constants 

We briefly recall the connection between the Lagrangian (jHJ) and the coupling con- 
stants gg^ Pl p 2 entering in Eq. (|T2*|) . After evaluating the traces, the Lagrangian 
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(JSJ) can be rewritten as a sum over all the decay channels S — > P1P2, where 
5 = {a^,a°, /c ± , /c , ko,cr,fo} and P = {7^, tt°, ?7, rf}: 

jCs^pp = CiSij (A^A'P) - c 2 S l3 (A'A'V 2 ) 

= (V2cx + ^=)(a° K-K + - K°Z°) + (ci + c 2 )(2a B n + n- + a B (ir ) 2 ) + .... 
v2 

= *s-> Pi p 2 SPiP2, (43) 

S,Pi,P 2 

where in the second line we wrote the a® to KK and the 03 to nn couplings explic- 
itly; they serve as illustrative examples of the adopted conventions. The coupling 
constants at this stage (without sum/average over isomultiplets) are simply given 
by g s ^p 1 p 2 = s ■ \s^p 1 p 2 , where s = 2/^/2 = v2 when P\ = P 2 , 1 otherwise. 

We have, for instance, g a o^ K+K - = (y/2ci + s/\j2c 2 ), 9 a o_^ K o = ~(V2ci + 



y/T/2c 2 ), ga B ^+n- = 2(ci + c 2 ), g aB ^ n o n o = \/2(ci + c 2 ) . Notice that by plugging 
these coupling constants into (|12|) we get the partial decay width for the correspond- 
ing channel, such as F a o^ K+K -. 

Let us then group the final states in corresponding isomultiplets: P = {n, K, 77, 77'}. 
One has to perform the sum over the members of the isomultiplets, for instance: 

dlo^KK = 9 2 al->K+K~ + 0fj>_S°jro ~> 9a°-+KK = a°-+K+ K~ (44) 

dlg^mv — QIb-^tt+it- + QIb-htW ~> 9a B ^mr = y^&TB-wr+TT- (45) 

The sign of the resulting coupling constant can be chosen arbitrarily for the 
isotriplet and isodoublets clq and k states. In the a[j to KK case we chose the sign 
of the charged decay channel K + K~ (which is opposite to the K K° one); Table 1 
is compiled following this convention. This however does not affect the decay rates 
expressed in (|12|) . Care is needed in the scalar-isoscalar sector because of mixing 
occurring in Eq. (|14|). such as g aB ^inr '■ the sign of g aB ^nn is taken to be the same 
as g aB _+ 7V+7r - and g aB ^^^- In the scalar-isoscalar sector the various contributions 
to the sum over final isomultiplets have the same overall sign, thus making the 
definition unambiguous. Table 2 and 3 follow this convention. 

As a last step an average over the initial isospin multiplets is performed, that is 
we consider S = {ao, k, a, fo}. For instance, Q 2 ao _^ K is given by: 

dlo^KK = 3 {/ll-tKK + 9 2 a +^KK + d^-iKlc) = ^l^KR- ( 46 ) 

The terms of the average are equal (this is true for each decay of ao and k\ for a 
and /o no average is needed). 
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Fig la 




Figure 1: Dominant (l.a) and subdominant (l.b) contributions to the transition 
amplitudes of a scalar tetraquark state into two pseudoscalar mesons. 
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